Electron transport in the four-lead two-impurity Kondo model: 
Nonequilibrium perturbation theory with almost degenerate levels 
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The eigenstates of an isolated nanostructure may get mixed by the coupling to external leads. 
This effect is the stronger, the smaller the level splitting on the dot and the larger the broadening 
induced by the coupling to the leads is. We describe how to calculate the nondiagonal density 
matrix of the nanostructure efficiently in the cotunneling regime. As an example we consider a 
system of two quantum dots in the Kondo regime, the two spins coupled by an antiferromagnetic 
exchange interaction and each dot tunnel-coupled to two leads. Calculating the nonequilibrium 
density matrix and the corresponding current, we demonstrate the importance of the off-diagonal 
terms in the presence of an applied magnetic field and a finite bias-voltage. 
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I. INTRODUCTION 

In equilibrium, the occupation numbers, or more gen- 
erally, the density matrix of the quantum states of an 
interacting electron system as a function of temperature 
T, chemical potential fi or magnetic field B, are given ex- 
plicitly in terms of the statistical operator. Out of equi- 
librium, even in a stationary situation involving e.g. a 
current flow through the system, the statistical operator 
is not known in general. Statistical expectation values of 
observables may however be calculated for noninteracting 
systems, as e.g. in the Landauer formula for the conduc- 
tance, or else in perturbation theory in the interaction 
using the method of nonequilibrium Green's functional. 

A particularly simple example of the effect of a finite 
current on the occupation of a quantum state is the po- 
larization of a spin-1/2 on a quantum dot in a magnetic 
field coupled by exchange interaction to the conduction 
electron spins in the leads. For sufficiently large chemical 
potential difference between source and drain electrodes, 
Ms — \ l d > 9^bB, where g^BB is the Zeeman splitting 
of the local spin levels, the occupation numbers are no 
longer determined by the thermal Boltzmann factors, but 
by a rate equation. More generally, it is a quantum Boltz- 
mann equation which describes the steady state of tran- 
sitions between the two local Zeeman levels induced by 
the available excess energy of electrons moving from the 
reservoir at higher chemical potential (source) to that 
with lower chemical potential (drain). These processes 
are mediated by the exchange interaction and the result- 
ing occupation numbers, or equivalently, the spin polar- 
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ization, can be very different from their thermal equi- 
librium values, even in the limit of vanishing exchange 
coupling. For example, the spin susceptibility, which in 
equilibrium obeys the Curie law \ ^ ^/T, is found to 
decrease with bias voltage V as x x 1/V at T <C (cf. 
Refs.SJl). 

In the case of a single quantum dot characterized by 
a spin-1/2, the density matrix of the local spin states is 
diagonal and the occupation numbers may be obtained, 
at least in lowest order in the coupling, by solving a rate 
equation. This is no longer the case in more complicated 
situations, when the density matrix is not even approx- 
imately diagonal in the basis of eigenstates of the iso- 
lated nanostructure or quantum impurity. The quantum 
Boltzmann equation therefore takes the form of a matrix 
integral-equation which can easily become numerically 
challenging. 

We show in this paper how the quantum Boltzmann 
equation may be solved approximately in a controlled 
way with the aid of a non-unitary transformation of 
the matrix Green's functions for the quantum impu- 
rity. This transformation serves to diagonalize the impu- 
rity (matrix-)spectral-function. By subsequently neglect- 
ing the broadening of these eigenstates of the cotunnel- 
coupled impurity, one again arrives at a simple rate- 
equation for the occupation numbers, valid to leading 
order in the cotunneling amplitude (e.g. second order in 
the exchange-coupling, i.e. fourth order in the tunnel- 
ing amplitudes) . This Bloch-Redfield^i 7 -^ type equation 
still involves off-diagonal terms, describing the transition 
amplitudes for the (voltage-)driven impurity system, but 
it remains far simpler than the full quantum Boltzmann 
equation and is readily solved numerically. The Bloch- 
Redfield equations for the reduced density-matrix of a 
quantum dot have been employed in the sequential tun- 
neling regime (cf. e.g. Refe.lsUlflD. Here we demonstrate 
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near the Fermi level, N(0), but with generally differ- 
ent chemical potentials \i n . The spin-operator s nm — 



FIG. 1: (color online) Double quantum dot setup: Two Kondo 
impurities, L, R, both represented by a spin-1/2 are coupled 
mutually by a spin exchange interaction K. The quantum 
dots are connected to two leads each, 1, 2 and 3, 4, respec- 
tively. 



how to establish similar equations in the cotunneling, 
or Kondo regime, starting from the quantum Boltzmann 
equation. 

As an example of a system where the off-diagonal en- 
tries of the density matrix become important, we study 
a system of two quantum dots in the Coulomb block- 
ade regime, each accommodating a single spin-1/2, mu- 
tually coupled by a spin exchange interaction, K . Each 
dot is contacted by a set of source and drain electrodes 
as illustrated in Fig. [TJ The eigenstates of the isolated 
two-impurity system are the singlet and triplet states of 
the two coupled spins. Nevertheless, cotunneling via the 
leads may mix the states, whereby the density matrix of 
the two-impurity system will acquire off-diagonal terms. 
As demonstrated below, this mixing can only occur in 
the presence of an applied magnetic field and for asym- 
metric couplings to the leads. The simpler case of zero 
magnetic field has already been analyzed in a previous 
publication^. For K = 0, the product states of the two 
spin-1/2 would be a natural basis, but for any finite K, 
perturbation theory should be performed with respect to 
the singlet /triplet basis. We will demonstrate how the 
singlet/triplet basis can be used, even for K = 0, as long 
as the off-diagonal parts of the density matrix are prop- 
erly taken into account. 



II. FOUR-LEAD TWO-IMPURITY KONDO 
MODEL 

We model this four-lead two-impurity Kondo model as 
illustrated in Fig. [TJ by the Hamiltonian 



(1) 



kncr 

+ KS L ■ S R - g^ B B ■ (Si + 9 R ) 

E ' s nm + z2 Jr™ S_R ■ S nm , 

n,m— 1,2 n,m— 3,4 

where n, m = 1,2,3,4 labels the leads, which are 
characterized by the same constant density of states 



k,k' 



C nk(r(V^'/2)Cr7 



denotes respectively the 



conduction electron spin (n = m) and the exchange tun- 
neling operator (n ^ to). Notice that we do not allow for 
charge transfer between the two dots. As we have demon- 
strated earlier—, the coupling to two pairs of source, and 
drain electrodes gives rise to a marked transconductance 
signal, reflecting the onset of cotunneling current through 
the left dot, say, when tuning the voltage over the right 
dot to match the exchange coupling, i.e. fi^ — ji^ = K . 

The eigenstates of the isolated two-impurity spin sys- 
tem are the singlet and triplet states 

\t+) = I T>l| T)/j, 

N = ^|(l T>L|l)^ + ll>L|T>fl), 

I*-) = I 1>h, 

\s) = ^ (I T>i|i>H-|i>i|T>fl). 

For K = 0, however, it would be more reasonable to use 
the product states 

\1) = \])l\[)r, \2) = \i) L \\) R , 

|3) = | T>i| T>«, |4) = | |)i| i) R , 

of the left and the right spin as a basis. Expressing the 
latter basis in terms of the former, the eigenstates with 
total spin quantum number S z = are seen to mix: 

|i) = IT)iU>« = ^(l*o> + |*», 

|2> = U>i|T)fl = ^(|*o)-|*)), (2) 

|3) = |< + ), |4> = |*->- 

In order to have a convenient representation of the 
spin operators in the singlet-triplet basis, we define a set 
of pseudo-boson (pb) operators {b^} = {s^ , t' + , tj, t_}, 
i.e. 7 € {s, i+, tg, £_}, to describe creation (annihilation) 
of a singlet state s^(s) or a triplet state rf,(t 7 )ii. The 
operators b^ span an infinite dimensional Fock space, 
which has to be projected onto the physical Hubert space, 
in which only single occupancy is allowed, i.e. Q — 
sKs + t\to + t+t + + fit- = 1- This constraint is enforced 
by adding a term XQ to the Hamiltonian and taking the 
limit A — > when calculating physical observables. 

The energy eigenvalues of the four states are 
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us = --K, wt = —K, 



and therefore 



K S L • S R = Ugsh 
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In terms of the pseudo bosons, the spin-1/2 operators at 
the left/right dot are given by 

S Z L/R = \{±sh Q ± 4s + t\t + - ttt.), (3) 



S 1/r = i S l/ R ) ] = l( ±sh ~ T A s + *+*o + 4*-)- ( 4 ) 
Or in compact notation, with a = L, R, 

S - - V fetT ,b , 

2 1 ' c*i77 7 ' 
7,7' 

where T Q;77 / is a vector of three 4x4 matrices, T x , T y , T z , 
defined by 



rpZ 

1 L/R — 



( ±1 

10 

±10 

V o o o -1/ 



1 L/R - 



±1\ 

Tl 1 

1 





where T+ = (T~)t = ^(T* + iTV). 

It is worth noting that the representation of Sl/r, 
Eqs. ((3]) and (j4]), does not include an s^s term, but 
only transition operators from a singlet to a triplet state. 
This implies, for example, that exchange-tunneling cur- 
rent cannot pass through either of the two quantum dots 
if the two-impurity system is in the singlet state. Inter- 
estingly, the excitation gap, K, to the current-carrying 
triplet states can be overcome by a finite bias across ei- 
ther of the two dots. This gives rise to a pronounced 
transconductance signal, which we have investigated ear- 
lier in Ref. 13. This work was restricted to zero magnetic 
field and hence avoided the problem of off-diagonal terms 
in the nonequilibrium correlation functions. 



III. DENSITY MATRIX AND OCCUPATION 
NUMBERS 

Assuming the exchange-tunneling to be weak, 
i.e. iV(0)J™ m <C 1, we shall determine the singlet-triplet 
occupation numbers by means of nonequilibrium pertur- 
bation theory. We employ contour-ordered Green's func- 
tions arranged in the matrix-form 



G 



G< + G r G> 
G< G> - G' 



satisfying the Dyson equation 



Gn G 



EG, 



where G and Go are the dressed and bare contour- 
ordered matrix Green's functions and £ is the self en- 
ergy. 

Whereas the information on the energy spectrum, in- 
cluding any shifts by external fields, is encoded in the re- 
tarded Green's function G r , the information on the ther- 
modynamic state of the system, i.e. the occupation of 



the energy levels, is contained in G < . Out of equilibrium, 
these two functions are not simply connected through 
a fluctuation-dissipation theorem and instead one must 
solve one more component of the matrix Dyson equation, 



Gn G K = S r G < 



£<G a . 



to 



An equivalent equation is obtained by applying G 
the second time argument in G (cf . Ref. ■ In contrast 
to the Dyson equation for G r , this quantum Boltzmann 
equation is a self-consistent equation for G < , since E < 
also depends on G < and since no unrenormalized Gq is 
known for an isolated quantum impurity. In the limit 
of vanishing coupling to the leads, G < is found to be 
independent of that coupling. Therefore the occupation 
number 



can take a finite limiting value dependent on the system 
parameters temperature, bias voltage and magnetic field 
in 0th order in the coupling to the leads (cf. also Refs. 0, 

i)- 

In the example of a two-impurity system, the pb 
Green's functions form a matrix in the singlet-triplet 
(ST) basis 7 = {s, t + , i , i.e. 



g r = 



(G r s 
1 
G 



G 



V 









r 

t+t+ 








fir 
u t t 









GLtJ 



where only the elements allowed by symmetry are shown. 

Notice that in the absence of any coupling to the leads, 
but assuming coupling to a heat bath, the quantum im- 
purity is in thermodynamic equilibrium and the Green's 
functions are given immediately as 



G<f\u) = -in B (w)A$(u), 
G>f )(cj) = -i[n B (u) + l]i4$(w), 

G 7 / a <^H 

4°7 ) (-) 



cj — uj-y — A ± iS ' 



2ttS(uj 



A), 



(5) 
(6) 

(7) 
(8) 



where lu 7 is the energy of the state 7 € {s, £+, to, i_} 
and it,b(oj) is the Bose distribution function. In the limit 
A — * 00 the Bose function at the position of the spectral 



peak, ui 



A, turns into a Boltzmann factor n 7 oc 



e -/3(^ T +A) ^ 1 The difference between the Bose and the 
Fermi statistics for the pseudo particles is then seen to 
vanish. Any term containing a product of two occupation 
numbers n^riy oc e~ 2/3X will be projected out at the end 
of the calculation. 



A. Pseudo-boson self energy 

The first order (in N(0)J™ m ) pseudo-boson self energy 
is proportional to the spin polarization of the conduction 
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electrons. It provides only a correction to the g-factor of 
the local spin and will be neglected in the following. 

The second order self energy, on the other hand, cor- 
responds to the diagram in Fig. [2] and reads 

S 7 ' 7 (ti,t 2 ) = --L Y a (T l7 T 2 ) (9) 

a=L,R 

Here the time variables t\ t% lie on the Schwinger con- 
tour, and a = L(R) for m,n — 1,2,(3,4). We have 
introduced the abbreviation, 

y Q (r 1 ,r 2 ) = 2 £ j™"Jr^?(n,r 2 ), (10) 

m,n=l,2,(3,4) 

with the summation variables depending on a and with 

^ % > k,k' 

being the conduction electron susceptibility. 

This expression is quite general and would be valid for 
any quantum impurity with internal states 7, where the 
matrices T l a would have to be defined accordingly. The 
self energy will have off-diagonal components (7 7^ 7/) 
if the conservation laws for spin allow so. Since a con- 
duction electron tunneling through one of the dots can 
either flip its spin or not, the accessible intermediate 
states may change the quantum number S z by 0, ±1. The 
second process has to flip the conduction electron spin 
back. Therefore, the self energy is diagonal in the quan- 
tum number S z . For the double dot considered here this 
leaves only one possible off-diagonal element, sto (and its 
hermitian conjugate), given by 

£ s i (n,T 2 ) = -7T (G tos (Ti,T 2 )[Y L (Tx > T2) +Y R (t 1 ,T 2 )] 

lb 

+ [G t _ i _(ri,r 2 ) - G t+ t + (Ti,T 2 )} 
x[Yl{ti,T2)-Y r (t 1 ,t 2 )]). (11) 

All other elements of the self energy are given in ap- 
pendix [5] We emphasize that the off-diagonal self en- 
ergy is finite only if two symmetries are broken simul- 
taneously: time reversal symmetry by a magnetic field 



7' / 13, 13' \ 7 



FIG. 2: Diagram for the second order pseudo boson self 
energy. Solid lines: conduction electron Green's functions, 
dashed lines: pseudo particles of the double quantum dot 
system. 



(t_l_ 7^ t—) and parity, i.e. the left-right symmetry (Yl 7^ 
Yr). 

This can be understood from the following simple ar- 
gument. Starting from the singlet state, a flipping of the 
left spin, say, causes a transition to a triplet state with 

S*\8)=±-j=\t T ), (12) 

which, upon a subsequent flipping of the left spin, makes 
a transition to either |s) or |io)- This is illustrated in 
Fig. [3] The two intermediate triplet states come with 
opposite signs, i.e. shifted in phase by 7r, and in the case 
of zero magnetic field these two alternative paths from 
\s) to I to) cancel and E sto vanishes. This is not the case 
for the diagonal component Y, ss in which the signs are 
squared and cause no cancellation. The observed depen- 
dence on the left-right symmetry, indicated in Fig. [3l 
arises in a similar way. 

Performing the analytical continuation to the real-time 
axis, one finds after Fourier-transformation in the relative 
time- variable, 

Y<(cj) =(-2tt) [(g 2 n + g 2 22 ) 2B(us) (13) 
+2. 9l2321 (B(lu + eV L ) + B(u - eV L ))] , 

Y<M =(-2tt) [(g 3 2 3 + .944) 2 B(w) (14) 
+2.934343 (B(uj + eV R ) + B{lu - eV R ))] , 

and <7 nm = N(0)J nm ,Y< ±R {u) = F<(w)±F<(w), where 
we have introduced the function 

B(x) = x ■ n B (x) = [coth(/3a;/2) - 1] x/2. 

From these correlation functions, the lesser component 
and the imaginary part of the retarded self energy are 



+/- 




-/+ 

FIG. 3: (color online) Illustration of the symmetry in the 
off-diagonal self energy £ s t - The two different paths from 
the singlet-state |s) to the triplet-state [to) over |t_) and |t+) 
come with opposite signs. Green(left) sign for interaction with 
left leads. Red(right) sign for interaction with right leads. 
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obtained as 



1_ f de_ 
16 7 2^ 

+ (G< f j e )-G< t+ ( e ))r<_ 7? (c- e )} 

(15) 

1 /" de 

16 7 2tt'' 

+ (A t _ t _(e)-A t+t+ (e))Y<_ R (e-u)}. 

(16) 



We observe that the off-diagonal spectral function obeys 
the sum rule 



r st0 M = MmE^ « -2-^ / ^l{A tos (e)y< +R (e-o;) 



did 

2^ 



A sio (o;) = ([ S ,4]) 



0. 



From the hermiticity condition it follows that 

[Gl t0 {u)Y = G a t0 M, 
and thus the spectral functions, 

A sto (uj) = A toS (uj), 
are identical. In addition one has 

[G^(t)]* = -G< a (-t) 

Re[G<»] = -Re[G<»] 
Im [G<»] =Im [G< S H] 



(17) 



(18) 



The diagonal elements of G^ 7/ are purely imaginary func- 
tions, and the off-diagonal elements share this property 
approximately 

Re[G<»] = -Re[G<»] «0 
^G<»«G<» (19) 

Thus the lesser Green's function is assumed to be sym- 
metric in analogy to the spectral function which was 
proven to be symmetric. It follows straightforwardly that 

E s t = ^tos- 



B. Retarded Green's function 

Assuming that the off-diagonal self energy is finite, we 
find by solving Dyson's equation Q r = [{Q r, ^)~ 1 — E r ] _1 , 



G 1 ' 



1 

dct 



/(GU)" 1 
det • G\ 






-'sto 





-'sto 








(GL)" 1 o 

det-GL t _/ 



EJL and det = 



where (G^ T ) = u> - uj 1 

(GU)- 1 ^.)- 1 "^) 2 . 

To lowest order in the coupling to the leads the re- 
tarded self energy is of 0(g 2 ), and the off-diagonal 
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FIG. 4: (color online) Spectral functions A ss ,A a t and At t 
versus the frequency uo/B for an exchange spin interaction 
K of the order of the level broadening, K — 0.05 ~ F sto . 
Further parameters are B = 1.0, §l = 0.1, qr = 0.2 and 
T — 0.001. The off-diagonal spectral function A sto has total 
spectral weight zero and thus contains negative values. 



Green's function is given by E£ t /det. Since the total 
spectral weight of the off-diagonal spectral function van- 
ishes, the function changes sign. It can be shown that it is 
approximately given by the difference of two Lorentzians, 
see illustration in Fig. 2J Since A sto (e) is not character- 
ized by a single peak, integration over e of a product of 
A sto (e) and a more slowly varying function B(e) may not 
be approximated as usual by taking B(e) at the position 
of the peak. 

In order to avoid having to deal with a not positive def- 
inite spectral function we may diagonalize the retarded 
Green's function matrix. The transformed matrix is 
given by 



\ Hl1^77 7'72 ^7l72> 

(Gil 








and 



V 



o 

t+t 






G 



22 



7 



G r t 

where G^ = l/(u> — and G£ 2 = l/(w — lu^) with 
1 



± y(cj t0 + ■S^JJ) -u s - ELM)' + 4 (S- t J ; 



This rotation is important only in the case when E£ t 
becomes of the same order of magnitude as the first term 
in the square root, which is proportional to the singlet- 
triplet splitting uj to ~lu s — K. The transformation matrix 
is given by 



U r 



(x\ -x\ 0\ 

10 

c r 2 x[ 

\ 1/ 



6 




FIG. 5: (color online) Spectral functions Au(uu) and A22(u) 
versus the frequency lu/B in the rotated space for the same 
parameter set as in Fig. |4] 



with the already normalized value of x\ = (x')i/J\f r and 



x 2 = {%')l/N r where Af r = \J {{x')^) 2 + {(x') r 2 f 



00; 1 «glh)- 1 -(g; (0 n)- 



and 



-e: + „. 



The transformation matrix U r has complex valued el- 
ements and its inverse is equal to its transpose. For 
K = 0, we find x\ — — x 2 — 1/V2 and the eigenstates 
of the system are given by the product states of the 
double quantum dot system, |1) = -^j(l^o) + I s )) an d 

|2> = ^(l*o)-|-». 

The transformation of the advanced Green's function 
is different, 

(Jjaylgajja _ ga ^ 

where the transformation matrix U a is defined in anal- 
ogy to the retarded Green's function with r replaced by 
a. The transformation matrix U a is related to U r by 
[jjajt = [zr]- 1 , so that [£ r ]t = G a . Since G^,G£ 2 are 
characterized by a single pole, the corresponding spectral 
functions show a single sharp peak and the commutation 
relations of the boson operators guarantee the integrated 
weight unity. In Fig. Owe show a typical example of the 
spectral functions Aa(u>) = — 21m G^, i = 1,2 . 



C. Lesser Green's function 

As mentioned above, the lesser Green's function has 
to be calculated self-consistently, e.g. in lowest order 
perturbation theory, thus in this model to second order. 



The lesser components of the Dyson equation are given 
by one of the equations 

(g*-)- 1 g< = i:<g a , 
g<(g a y 1 = g r z<, 

where we neglected the boundary terms. Subtracting 
these two equations, G < is found to obey the quantum 
Boltzmann equation 2 , 



(ST 1 Q K - Q K [Q'T 1 = ^ < g a -g r v 



< r.o. n.r v< 



(20) 



If we neglect the off-diagonal terms in Eq. (f2"0"|) we have 
to solve the following equations, 



(G r ) - {G a y 



77 



G< = S< 7 [G a ~ G r l n 



G< 



77 



r 77 (w) 



which is a self-consistency equation for the occupation 
number n~ ( — i S 77 (u; 7 )/r 77 (w 7 ). Since the broaden- 
ing of the spectral function j4 77 is the smallest energy 
scale in the problem, we neglect the frequency depen- 
dence and consider only the on-shell occupation num- 
bers. In this lowest order approximation the quantum 
Boltzmann equation (|20|) is essentially a rate equation 
(see appendix [Bl for explicit expressions of the occupation 
numbers), albeit one which now includes the off-diagonal 
components of the density matrix. 

We will now demonstrate how this approach fails in the 
example of the magnetization. For K = 0, the left and 
the right spins as defined in Eq. (J3]) are good quantum 
numbers and we can define the corresponding magneti- 
zations 



M L/R = 2{Sl /R ) = (t\t + 



under the assumption that (sHo + t' Q s) = 0. This leads 
to 



M 



L/R 



Y<(B) + Y<(B)-Y<(-B) 



-B) 



Y< (B) + Y< (B) +Y<(-B) + Y< (-B) ' 



from which the magnetization on the left quantum dot is 
found to depend on the magnetization of the right quan- 
tum dot and vice versa, even though the dots are com- 
pletely decoupled. For a single Kondo impurity, however, 
the correct result is£^ 



Mr 



Y<(B)-Y<(-B) 
Y<(B)+Y<(-B) 



(21) 



The off-diagonal components are not important in the 
case of left-right symmetry and indeed for Yl = Yr we 
obtain the correct expression for Ml,. Also, at B = and 
consequently Mj, = the two different results coincide. 
The difference of the two above results for Ml can be 
traced back to the unjustified neglect of the off-diagonal 
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average (sHq + t^s). We show in the appendix [Cl that 
including the effect of off-diagonal terms, by employing 
the transformation defined above, we recover the correct 
result for the magnetization Ml, Eq. (f2~Tj) . From this 
calculation it is obvious that the off-diagonal terms con- 
tribute to the same order as the diagonal terms. 

In the case of finite exchange interaction K the calcu- 
lation including the off-diagonal contributions becomes 
cumbersome, since one has to solve a self-consistent sys- 
tem of integral equations. The solution can be much 
simplified in the transformed basis introduced above. 

To rotate the Quantum Boltzmann equation we mul- 
tiply Eq. (20]) with U r from the left and [U" 1 )- 1 = (J7 r ) t 
from the right. Thus we find the transformed quantum 
Boltzmann equation, 



g< (g a 



s< g a - g r sr 



(22) 



where t< = C/ r E < ([/ Q )- 1 = £/ r £ < (L/ r ) t and §< = 
U r g<(U a )- 1 = U r g<{U r Y. After the transformation 
all retarded and advanced Green's functions are diago- 
nal and we obtain a single equation for every entry of 
the lesser Green's function. There is still a finite off- 
diagonal element 2 ^ in the rotated basis, and the big 
advantage over the initial formulation is the fact that the 
spectral functions appearing on the r.h.s. of Eq. (j2"2"|) are 
now all positive definite and may be approximated by 
delta functions of weight unity. As for the corresponding 
real parts of g a,r we adopt the usual assumption that 
after frequency integration they may be neglected. From 
Eq. (|22| Q < is then found as a sum of delta functions 

— 7 

with weights N to be determined self-consistently, 



-i N'2tt5{lu-uj^). 

7=1, 2, t + ,t- 



This expression, when substituted into the lesser self en- 
ergy, leads to a linear combination of weight factors. The 
quantum Boltzmann equation reduces to a set of linear 

~7 

homogeneous equations for N , which, together with the 
normalization condition, s^s + t^to + t\t+ + t_t- = 1, 

~7 

may be solved to give the weight factors. Note that N 
are matrices in the local Hilbert space defined by states 
(l,t+,2,i_) with nonzero components N^ +t+ , N t 
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t_t_ 



JVJi, N? 2 and JV£. 
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The expectation value (sHq + i^s) is obtained from 



st 



as 



(sh + tls) = i \ — [G<» + G<»] 



2tt 

du> r 

2^ 



i) (pfi,i) 



(2,2) 



+ {XyX'i X2X2) f^(l 
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+ (x{xl - x$x a 2 ) (7V7 2 - Nl) 
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FIG. 6: (color online) The expectation value of (sHq +t s), 
i.e. the contribution of the off-diagonal Green's function is of 
the order of the diagonal contributions as illustrated here for 
the cases of K/F 3 t — 1, 0.2, 0. Further parameters of the plot 
are B = 1.0, g L = 0.1, g R = 0.2 and T = 0.001. 



A numerical evaluation of (sHq + ijs) is shown in Fig. [5] 
For the parameters chosen, it is seen to be of the same 
order of magnitude as for example the non-equilibrium 
magnetization, and is thus comparable to the diagonal 
occupation numbers. In Fig. we compare the results for 
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FIG. 7: (color online) The magnetization Ml of the left quan- 
tum dot is strongly influenced by a voltage applied over the 
left dot. For K ~ Y s t , the magnetization, Mr, over the 
right dot shows only minor deviations from the thermody- 
namic value. Neglecting the off-diagonal contributions would 
give Ml = Mr = M to t/2. Further parameters of the plot are 
B = 1.0, g L = 0.1, g R = 0.2 and T = 0.001. 

the magnetization with, and without off-diagonal compo- 
nents, in the parameter regime K ~ T, where the effect 
of the off-diagonal components was argued to become im- 
portant. The magnetization of the left/right spin is given 
by the sum/difference of the total magnetization and the 
off-diagonal contributions 

M L/R = 2(S L/R ) = ((S L + S R ) ± (S L - S R )), 
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where (S L + Sr) = n t+ - n t _ = M tot /2 and (S L - S R ) = 

(sHq + i^s). For K = the magnetization of the right 
quantum dot should not depend on the voltage applied 
to the left quantum dot. Therefore the off-diagonal ex- 
pectation value, (sHq +£gs), has to compensate 100% of 
the voltage-dependent part of Mtot/2- For the parameter 
regime in Fig. [7J K ~ V st , the compensation from the 
off-diagonal contribution is already of the order of 75%. 



IV. CALCULATION OF THE 
NON-EQUILIBRIUM CURRENT 

The non- linear conductance dl/dV is governed by the 
voltage-dependence of the occupation of states. For in- 
creasing voltage V, the conductance will have a step 
when the energy supplied by V allows the occupation 
of an excited state. The voltage-dependence of the level 
occupations will change the step to a cusp, and higher 
order in the perturbation series will change this to a log- 
arithmic nonequilibrium Kondo peak, cut-off by the spin- 
dependent relaxation rate. 

To second order in the exchange-tunnel coupling the 
current through the left quantum dot is given by ^ 14 i 15 i 16 



7r e 
8h 



#12321 



eV L )-B{u + eV L )\, 



where X^q D (u>) is the susceptibility of the double quan- 
tum dot system 



x dqd( w ) 



all 



*Tr 
2tt 



G<{e)TlG>{e + u)Tl 



This expression has to be modified according to the ro- 
tation in the basis states. Therefore we use the rep- 
resentation Q y = Q r — Q a + Q < , together with the 
fact that terms containing more than one factor of Q < 
are projected out. The lesser Green's function is given 
by Q < — U a G < (U r )~ 1 and the retarded and advanced 
Green's functions, Q r and Q a , have to be transformed ac- 
cording to G r / a = (U r / a )~ 1 g r / a U r / a . The two limiting 
cases E s t = and K = are discussed in appendix [D] 
The off-diagonal contributions show a significant effect 
also in an intermediate regime. 

In Fig. [5] the differential conductance is plotted for the 
parameters K/T sto ~ 1,0.2,0, and for K ~ T sto both 
with and without the correction caused by a finite expec- 
tation value (sHo+t s). One observes a significant differ- 
ence, in particular near threshold. Since the current ex- 
pression depends sensitively on the non-equilibrium occu- 
pation numbers, the physically relevant results are given 
only for the correct occupation numbers. If the finite con- 
tribution of (sHq + 4 s ) is neglected, the current through 
the left quantum dot depends on the total magnetization 
of the double quantum dot system, such that e.g. a fi- 
nite voltage on the right quantum dot would affect the 
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FIG. 8: (color online) Differential conductance (in units of 
(7r/8) (e 2 /h) 312321) in the case of small exchange interaction 
K/Y a t — 1, 0.2, 0. For K ~ r st() the results including or ne- 
glecting the off-diagonal contributions are shown. The further 
parameters are chosen identically to the previous figures. 



current through the left, although the two quantum dots 
are decoupled when K = 0, see also discussion in the 
appendix [D] 



V. CONCLUSION 

The state of an isolated nanostructure is determined 
by the specification of the occupation of the eigenstates 
of the system. In other words, the density matrix (the 
lesser Green's function integrated over frequency) of the 
isolated system is diagonal in the basis of eigenstates. 
Coupling of the nanostructure to reservoirs will in gen- 
eral lead to a change of the density matrix. This change 
may involve the appearance of off-diagonal elements in 
the density matrix. These off-diagonal elements gener- 
ically have a more complex frequency dependence than 
the diagonal terms. While the diagonal terms of the den- 
sity matrix have a spectral function characterized by a 
single narrow peak, and a spectral weight to be inter- 
preted as the occupation number of the state in question, 
the spectral functions of the off-diagonal elements of the 
density matrix have positive and negative parts and to- 
tal spectral weight zero. Nonetheless these off-diagonal 
terms may be as important as the diagonal ones, as we 
demonstrate in the example of a double quantum dot 
system in a magnetic field. 

We have presented a systematic method of how to deal 
with this problem, by introducing the two sets of eigen- 
states of the retarded and advanced Green's functions. 
In terms of these eigenstates the quantum Boltzmann 
equation may be solved in the usual approximation of as- 
suming the spectral functions to be delta functions. The 
method is generally applicable, but it is demonstrated 
here in the example of a minimal model, where the prob- 
lem of off-diagonal elements of the density matrix arises: 
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a double quantum dot system coupled by spin exchange 
interaction in a magnetic field. 

In this case, the eigenstates of the isolated double-dot 
system are the singlet and triplet states, even for arbi- 
trarily small exchange interaction K. However, when the 
leads are coupled to the dot system, and K is of order, 
or less than the level broadening on the dot system in- 
duced by the leads, the eigenstates of the coupled system 
- quantum dots plus leads - approach the product states 
of the spins 1/2 of the individual dots. The transition 
in the character of states as the exchange coupling K is 
varied is captured perfectly by the representation in the 
rotated basis proposed here. Our method thus allows to 
avoid the time-consuming numerical solution of the full 
frequency dependent quantum Boltzmann equation. 

The method is quite general and can be applied to a 
wide range of quantum-impurity problems. This is par- 
ticularly relevant in multi-orbital problems such as car- 
bon nanotube quantum dots and single-molecule tran- 
sistors involving smaller conjugated molecules, possibly 

I 



£ ss (ti,t 2 ) 

S toto (Tl,T 2 ) 
S t±t± (Tl,T 2 ) 



and 

S s t (ri,T 2 ) 

St 0S (n,T2) 



APPENDIX B: SOLUTION OF THE QUANTUM 
BOLTZMANN EQUATION WITHOUT 
OFF-DIAGONALS 

If we neglect the off-diagonal elements in Eq. (|2H)) we 
have to solve the following equations, 

which corresponds to a self-consistency equation for the 
occupation number n 7 = i £ 77 (w 7 )/r 77 ((jJ 7 ), when ne- 
glecting the frequency dependence and considering only 
the on-shell occupation numbers. All other quantities 
will change on a larger energy scale than the broaden- 
ing of the spectral function A 7 , therefore we are allowed 



acting as high-spin impurities. 
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APPENDIX A: DERIVATION OF THE SELF 
ENERGY 

From Eq. © we find the following self energies (a com- 
mon prefactor of -1/16 is implied), 



(Al) 
(A2) 

(A3) 



(A4) 
(A5) 

I 

to approximate it by a (5-function. This simplifies the 
frequency integration in the self energies, leading to a 
set of homogeneous linear equations for n 7 . These equa- 
tions are closed by imposing the normalization condition 
^n 7 = 1. The solution is 

7 

AA 7 

where 

Zt = Y< +R (-K) + Y< +R (B) + Y< +R (-B), 
Zt ± = Y< +R (-K ±B)+ Y< +R (±B), 

and 



: (G Mo (ti,t 2 ) + G t+t+ (ri,r 2 ) + G t _ t _ (n, r 2 )) Y l+r (t 1 ,t 2 ), 

- (G ss {t 1i t 2 ) + Gt + t + (ri,r 2 ) + G t _ t _ (n, r 2 )) Y l+r (t 1 ,t 2 ), 

- (G ss {ti,t 2 ) + G to t (Ti,T 2 ) + G t±f± (ri,T 2 )) Y l+r (tx,t 2 ) 

T (G s t (Tl,T 2 ) + G ta s(Tl,T 2 ))Y L - R (Tl,T 2 ), 



(G f _ t _(ri,r 2 ) - G t+t+ (ri,r 2 )) Y l ^ r {ti,t 2 ) + G toa (n, t 2 )Y l +r(ti, t 2 ), 
(G t _ t _(T 1 ,T 2 ) - G t+t+ (r 1 ,r 2 )) Y l ^ r (tx,t 2 ) + G sto (n, t 2 )Y l+r (t 1 , t 2 ). 
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Af s =Z t0 Z t+ Z t _ - (Zt_ + Z t+ ) Y< +R (B)Y< +R (-B), 

Mo =Z t+ Z t _Y< +R (K) + Z t+ Y< +R (-B)Y< +R (K + B) + Z t _Y< +R (B)Y< +R (K - B), 
Af t± =Z tT Z tQ Y< +R (K T B) + Z tT Y< +R ( T B)Y< +R (K) 

+ Y£+r(tB)Y< +r ( T B)Y< +r (K ±B)- Y< +r ( T B)Y< +r (±B)Y< +r (K t B). 



APPENDIX C: QUANTUM BOLTZMANN 
EQUATION IN THE LIMITING CASE K = 

In the special case of K = the ss and toto components 
have the same energy, u s = u>t = 0. Their zeroth order 
contributions are identical and they will be the same in 
every following order of the calculation. The same argu- 
ment holds for sto and tos. 

Thus we have to solve only the two equations following 
from Eq. (3DJ> 



is of the same order as the diagonal contributions like 
(t_t- — tj_£+). Please note, that within this calculation 
the occupation numbers, ni = n^nn^ are given as prod- 
uct states. Although the two quantum dots are decoupled 
in the case of K = 0, the solution for the occupation num- 
bers of the product states contains information about the 
left and right quantum dot simultaneously. 
Finally, we find for the magnetization 
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^0*0 



'-'sto ^t o t - 



tos) 



G? 



+ ((s + t ) t (s + to)- (to-8)Hto 



0) 



=n L] n m - n Ll n Rl 



Y L <(B) 



L 



Yr<(- 



B) 



Y L <(B) 



The sum and difference of these two equations, using 
G s t = Gt s and G ss = Gt t m the special case of K = 0, 
are 



(G r ss ) 1 ± (G r sto ) - (G" s ) 1 =f (Gg tg ) 
= [^f s =t ^ft ] [Gss ± G°t — GJ S =f G^ t(J ] 



[G 



For zero exchange interaction ii" = we have to do per- 
turbation theory in the rotated eigenspace 



\s) 
to) 



where the new states |1) and \2) correspond to product 
states and are orthonormal. Again using the strongly 
peaked nature of the spectral functions we find a set 
of linear equations for the occupation numbers, which, 
amended with the normalization condition, give the so- 
lution 



n 2 



Y<(-B)Y<(B) 



(Y<{B)+Y<{-B)) (Y<(B)+Y<(-B)Y 
Y<(B)Y<(-B) 



(Y£(B) 



Y L < 



n t . 



(-B))(Y<{B)+Y<{-B)y 

Y<( T B)Y<( T B) 

(Y<(B) +Y<{-B)) (Y<(B)+Y<(-B)Y 



The discussion in the main text shows, that the contri- 
bution of the off-diagonal elements 

(sho + t' s) = n% - n 2 



APPENDIX D: EXPLICIT CURRENT 
EXPRESSIONS 

The general expression for the current will not be given 
here, but we like to discuss briefly two limiting cases. If 
the off-diagonal contributions are zero, for example for 
B = or left-right symmetry, the current is given by 



(^) 2 {\g 2 i2 [3eVb + (n. 



+ [Ji s - 
+ (n s - 
+ ("-to 



Tl t _ 

- n t _ 



where the function 
1 



-n t0 )F 3 (K,V L ) 

F 3 (K + B,V L ) 

F 3 (K~B,V L ) 

-n t0 +n t+ )F 3 (B,V L )] , 
(Dl) 



1 



F 3 (x, V) = -{x - eV) coth -/3(x - eV) 



[x + eV) coth -f3(x + eV) 



is asymmetric in the voltage V and also in the energy x. 
The function F 3 (x, V) describes the non-linear behavior 
of the differential conductance. 

In contrast, for the case of K = 0, the current is given 

by 



[3eV L 



(n t+ - n t _ 



ni-n 2 )F 3 (B,V L ) 
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This result is identical to the current for a single quantum 
dot in non-equilibrium^, since n t+ — n t _ + n\ — n 2 = M^. 
The current depends only on the magnetization of the 
left quantum dot. Taking the limit of K — > in the 
expression (ID1[) . the current is proportional to the total 



magnetization (n t+ —n t _) and thus dependent on proper- 
ties of the right quantum dot, which is obviously wrong. 
Please compare with the discussion of the magnetization 
in section UlI CI 
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